Atomic quantum dots coupled to BEC reservoirs 
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We study the dynamics of an atomic quantum dot, i.e. a single atom in a tight optical trap 
which is coupled to a superfluid reservoir via laser transitions. Quantum interference between the 
collisional interactions and the laser induced coupling to the phase fluctuations of the condensate 
results in a tunable coupling of the dot to a dissipative phonon bath, allowing an essentially complete 
decoupling from the environment. Quantum dots embedded in a ID Luttinger liquid of cold bosonic 
atoms realize a spin-Boson model with ohmic coupling, which exhibits a dissipative phase transition 
and allows to directly measure atomic Luttinger parameters. 



A focused laser beam superimposed to a trap holding 
an atomic Bose-Einstcin condensate (BEC) l] allows the 
formation of an atomic Quantum dot (AQD), i.e., a single 
atom in a tight trap 0, Q which is coupled to a reservoir 
of Bose- condensed atoms via laser transitions. This con- 
figuration can be created, e.g., by spin-dependent optical 
potentials Q, where atoms in the dot and the reservoir 
correspond to different internal atomic states connected 
by Raman transitions. Atoms loaded in the AQD will 
repel each other due to collisional interactions. In the 
limit of strong repulsion, a collisional blockade regime 
can be realized where either one or no atom occupies the 
dot, while higher occupations are excluded. Below we 
will study the dynamics of such an AQD coupled to a 
BEC reservoir: as the key feature we will identify the 
competition between two types of interactions, namely 
the coupling of the atom in the dot to the BEC den- 
sity fluctuations via collisions, and the laser induced cou- 
pling to the fluctuating condensate phase. Depending on 
the choice of interaction parameters, they can interfere 
destructively or constructively, providing a tunable cou- 
pling of the dot to the phonons in the condensate in the 
form of a spin-Boson model [jj. In particular, an essen- 
tially complete decoupling of the dot from the dissipative 
environment can be achieved, realizing a perfectly coher- 
ent two-level system. This interference and tunability of 
the coupling of the dot to the environment occurs for 
condensates in any dimensions. A particularly interest- 
ing case is provided by a ID superfluid reservoir 0, i.e., 
a bosonic Luttinger liquid of cold atoms pj, where the 
system maps to a spin-Boson model with ohmic coupling. 
The tunable dot-phonon coupling then allows the cross- 
ing of a dissipative quantum phase transition 0, 0] , and 
can serve also as a novel spectroscopic tool to measure 
directly atomic Luttinger parameters. 

Let us consider cold bosonic atoms with two (hyper- 
fine) ground states a and b (Fig. . Atoms in state a 
form a reservoir of atoms in a superfluid phase, held in a 
shallow trapping potential V a (x) . The AQD is formed by 
trapping atoms in state & in a tightly confining potential 
Vb(x) produced, e.g., by a focused laser beam induced 
potential or by a deep optical lattice potential which is 



FIG. 1: Schematic setup of an atomic quantum dot coupled 
to a superfluid atomic reservoir. The Bose-liquid of atoms in 
state a is confined in a shallow trap Va(x). The atom in state 
b is localized in tightly confining potential 14 (x). Atoms in 
state a and b are coupled via a Raman transition with effective 
Rabi frequency Q. A large onsite interaction Ubt > allows 
only a single atom in the dot. 



only seen by atoms in state b Within the standard 
pseudopotential description, the collisional interaction of 
atoms in the two internal levels a,f3 = a,b is described 
by a set of coupling parameters g a p = Airaaph 2 /m with 
scattering lengths a a p and atomic mass to. We as- 
sume that the reservoir atoms are coupled via a Raman 
transition to the lowest vibrational state in the AQD, 
where spontaneous emission is suppressed by a large de- 
tuning from the excited electronic states. Thus, follow- 
ing arguments analogous to those in the derivation of the 
Bose-Hubbard model of cold atoms in an optical lattice 
0, we obtain an effective Hamiltonian, 

H b + H ab = (-HSo + g ab J dx |V fc (x)| 2 p a (x)^ S+S (1) 

+ -^-Wbb + J dx^(* a (x)V'6(x)6 t +h.c), 

Here 'l'a(x) is the annihilation operator for an atom a 
at the point x, while /5 a (x) = #t(x)>I' a (x) is the asso- 
ciated operator for the density. The operator b destroys 
a &-atom in the dot in the lowest vibrational state with 
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wave function ^(x). The first term in Eq. includes 
the Raman detuning So, and the collisional interactions 
between the &-atoms with the reservoir. The second term 
describes the onsite repulsion Ubb ~ gbb/if > between 
atoms in the dot with 4 the size of the ground state 
wave function. The last term in Q is the laser induced 
coupling between a and b atoms with effective Rabifrc- 
quency fl. In writing JJJ we exclude coupling to higher 
vibrational states in the dot, assuming that these states 
are offresonant Q. 

At sufficiently low temperatures the reservoir atoms in 
a form a superfluid Bose liquid with an equilibrium liquid 
density p a . The only available excitations at low energies 
are then phonons with linear dispersion uj = t> a |q| and 
sound velocity v s . With the assumption that the num- 
ber of condensate atoms inside the dot is much larger 
than one, n a — p a l% >> 1, i.e. lb is much larger than the 
average interparticle spacing in the BEC reservoir, the 
quantum dot is coupled to a coherent matter wave and 
the Bose-field operator can be split into magnitude and 
phase, ^(x) ~ p a (x) l ^ 2 e~^ x K The dimensionless pro- 
portionality coefficient is the bare condensate fraction, 
which depends on non-universal short distance proper- 
ties of the Bose-liquid. This representation does not 
require the existence of a true condensate and can be 
also used to describe both 2D and ID superfluid systems 
[ill ]. The density operator can be expressed in terms of 
the density fluctuation operator IT: p a (x) = p a + II(x), 
which is canonically conjugate to the superfluid phase 
<f>. In the long wavelength approximation the dynamics 
of the superfluid is described by a (quantum) hydrody- 
namic Hamiltonian |llj ] 
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</x[ — /> s |v«£| 2 (x) + ^n 2 (x) 
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where p s is the density of the superfluid fraction (at zero 
temperature p s — p a ). The quadratic Hamiltonian (J2J 
is easily diagonalized by introducing standard phonon 
operators 6 q via the following transformation: 
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(6q - bU) 
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with V the sample volume. Accordingly, the Hamilto- 
nian J2J takes the form of a collection of harmonic sound 
modes: H a = hv s J2q |q|^q^q- Since the excitations of 
a weakly interacting Bose-liquid are phonon-like only for 
wavelengths larger than the healing length £ > lb, the 
summation over the phonon modes is cutoff at a fre- 
quency LOc = U s /£ ~ gaaPa/h- 

Iii the following we consider the collisional blockade 
limit of large onsite interaction Ubb, where only states 



with occupation rib = and 1 in the dot participate in 
the dynamics, while higher occupations are suppressed by 
the large collisional shift. This situation and its descrip- 
tion is analogous to the Mott insulator limit in optical 
lattices 0, 0, 0] . The requirements are that the Ra- 
man detuning and Rabi couplings are much smaller than 
Ubb/h, which provides a small parameter to eliminate the 
states with higher occupation numbers perturbatively. 
As discussed below, a Feshbach resonance can help in 
achieving this large Ubb limit 0|. Thus, the quantum 
state of a dot is described by a pseudo-spin-I/2, with the 
spin-up or spin-down state corresponding to occupation 
by a single or no atom in the dot. Using standard Pauli 
matrix notation, the dot occupation operator wb is then 
replaced by (1 + cr z )/2 while b< — > c + . Furthermore, the 
dominant coupling between the AQD and the superfluid 
reservoir arises from the long wavelength phonons. For 
wavevectors |q|4 -C 1, the phonon field operators may be 
replaced by their values at x = 0. Neglecting the density 
fluctuations in the Raman coupling (see below) and an 
irrelevant constant, the Hamiltonian Q is simplified to 
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9ab 



hA 



2 H(0)j a z + — (a + e-^+h.c.). 

(5) 

Here A ~ fln a ^ 2 is an effective Rabi frequency with 

1 /2 

n a the bosonic enhancement factor due to the BEC 
reservoir, and with a proportionality coefficient which 
depends on the bare condensate fraction and the ex- 
plicit form of the wavcfunction -i/v The form of the 
Rabi coupling in JSJ only applies for A <C u> c , which 
is the interesting regime in the spin-Boson model dis- 
cussed below. In (|5jl the detuning has been renormalized 
to include a mean field shift and a shift due to the vir- 
tual admixture of the double occupied state in the dot, 
-M + g ab pa + (hA) 2 /(-2U bb ) = -hS. 

Moreover, the validity of the above model requires 
a strong collisional interaction of atoms in the AQD, 
.966 3> g a a- This follows from the inequalities n a 3> I, 
A -C lo c and the single occupancy condition HA <C Ubb- 
A magnetic or (Raman laser induced) optical Feshbach 
resonance in the b channel will assist in achieving this 
limit 13]. These resonances arise from coupling to a 
bound molecular state in an energetically closed colli- 
sional channel. In the case of an optical Feshbach res- 
onance, for example, the onsite interaction due to the 
laser induced Raman coupling of two b atoms in the dot 
to a molecular bound state leads in Eq. to the re- 
placement Ubb —> ^6& eS ' ) = Ubb + g 2 /S m , where the second 
term describes the resonant enhancement with g an ef- 
fective Raman Rabi frequency and S m the detuning from 
the molecular resonance is (valid for g < \S m \). A finite 
lifetime of the molecular state, e.g. due to collisions with 
a atoms, introduces a width S m — > S m — iT/2 |2l|. Thus 



for detunings \S m \ T we have Ubb — * U { 
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with 7 m = (g 2 /S 2 n )T < U b % cs) . Returning to the Hamil- 
tonian (JjjJ we see that besides the resonantly enhanced 
onsite interaction we have a non-Hermitian loss term 
~ (hA/U bb s ') 2 'Y m which is strongly suppressed in the 
collisional blockade limit. 

Eventually after a unitary transformation H = 
S^iHa + Ht + Ha^S with S = exp(-cr z ^(0)/2) the dy- 
namics of the AQD coupled to the phonons of the super- 
fluid reservoir is described by a spin-Boson type Hamil- 
tonian [j| 



H 



HA 



q 




The coherence in the reservoir is crucial for the most 
important feature of our cold atom version of the spin- 
Boson model (©, namely the fact that the collisional in- 
teractions and those arising from the coupling of the Rabi 
term to the condensate phase add coherently in the am- 
plitudes of the total phonon coupling 
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In particular, we see that for a repulsive inter-species in- 
teraction, g a b > 0, both contributions interfere destruc- 
tively: the effect of the phonon excitation in a laser- 
driven transition a <-> b can be precisely cancelled by 
the change in the direct "elastic" interaction between the 
liquid and the b atoms, as described by the first term in 
Eq. ifHjl. In a weakly interacting gas mv 2 — p a g a a and 
thus the coupling constants A q vanish at g a b = g aa . At 
this special point, we have thus formed the analogue of 
a "charge" qubit of solid state physics ^ij, i.e., a coher- 
ent superposition of occupation and non-occupation of 
the dot, which - at vanishing detuning - exhibits perfect 
Rabi oscillations of the AQD's occupancy. Notice that 
the simple form Q is valid only for |q|Z& <C 1, while 
for larger values of the momentum the interaction starts 
to decrease proportional to the Fourier transform of the 
wavefunction tpb- 

We discuss now the importance of the neglected den- 
sity fluctuations in the Rabi term and their effect at 
the decoupling point. After the unitary transforma- 
tion above, these fluctuations give rise to a perturba- 



tion W 



ha n(o) 

2 2p a 



a x which is small compared to the col- 



lisional interaction in Eq.JSJ provided that HA -C p a gab- 
Now unless g a b is much smaller than g aa , this condition 
is precisely equivalent to the condition A<w c discussed 
above. In particular, if we are at the decoupling point, 
the only remaining interaction with the phonon bath is 
given by H' . For 6 = the Hamiltonian can then be 
rewritten as 



H = — 



HA 



1 - 



n(Q) 

2 Pa 



(8) 



This is an independent Boson model which can be di- 
agonalized exactly [lij . Since the phonons now no 
longer couple to the b atom occupation a z and, more- 
over, the coupling constants are proportional to the small 
parameter A/u c , one obtains perfect Rabi oscillations 
except for a small reduction in amplitude by a factor 
cxp — 7(A/cj c ) 2 ps 1, where 7 < {paO^a) 1 ^ 2 ^ 1 is pro- 
portional to the small gas parameter. 

Let us now turn to discuss the properties of system 
when it can be described by the Hamiltonian Eq. (jSJ. 
The system is characterized by the effective density of 
states 



J(w) 



J2 ^(uj - Wq) = 2aLU S 



(9) 



where a ~ (gabPa/mv^ — l) 2 is the dissipation strength 
due to the spin-phonon coupling and D = s the dimen- 
sion of the superfluid reservoir. In the standard termi- 
nology 5( , s = 1 and s > 1 correspond to the ohmic and 
superohmic cases, respectively. In the superohmic case, 
the resulting dynamics of the AQD is a damped oscil- 
lation at vanishing detuning 5 = 0, consistent with the 
result of simple Bloch equation analysis. It may be ob- 
served by following the population of the atoms b in the 
presence of the laser coupling. The associated frequency 
A and damping T can also be obtained by measuring a 
weak field absorption spectrum, which would exhibit the 
oscillations with frequency A as a line splitting, and T as 
the linewidth. 

A much richer dynamics appears for ohmic dissipa- 
tion (Fig. EJ. In this case, the system exhibits a zero 
temperature dissipative phase transition, as a function 
of the dissipation strength, at a critical value a c = 1 for 
A <C 0J C H- I n the symmetry broken regime a > a c the 
occupation probability of the b atom will exhibit a finite 
jump from (1 — m s )/2 to (1 + m s )/2 as the detuning 5 
is changed across zero. The spontaneous polarization m. 
is a function of a approaching m s ~ 0.9 for a — > |l8| 
and m s = 1 for a ^> a c . As a result the b occupation 
probability is almost unity for any a > a c . Instead in the 
regime a < a c , at vanishing detuning, the average pop- 
ulation of the b atoms is 1/2. In particular, for a < 1/2 
one has damped Rabi oscillations. In terms of the char- 
acteristic frequency scale A r (a) = A(A/uj c ) a ^ 1 ~ a \ the 
effective Rabi oscillation frequency A and damping rate 
r are given by A = cos 77 • A r (a) and T = sin 77 • A r (a) 
with 77 = 7ra/(2(l — a)) p. This result holds as long as 
T < T* = KA r /a. At higher temperatures the dynamics 
is incoherent and no oscillations should be visible. For 
1/2 < a < 1 the Rabi oscillations disappear and the 
behaviour is completely incoherent |l6|. Only numerical 
results are available for the dynamics in this regime 
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The ohmic spin-Boson model is achieved by embed- 
ding the AQD in an atomic quantum wire as realized 
recently in |19|. For transverse harmonic trapping both 
the temperature T and the chemical potential need to be 
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smaller than the frequency of the transverse confinement 
u>±. The proper description of a ID superfluid in terms 
of the hydrodynamics Hamiltonian (J2J is based on the 
Haldane-Luttinger approach 7]. The definitions of the 
phase and the density fluctuations in terms of the phonon 
operators are given by Eas. 1)3141) with v s — tmp a /(mK), 
where p a ~ p a l± is the ID density of the cloud and 
li_ = {h/muj^)^ 1 / 2 the transverse ground state size. The 
Luttinger parameter K depends on the density of the 
liquid through the combination 7 a = mg aa /h 2 p a - The 
corresponding ID interaction constant g aa is related to 
the 3D scattering length a aa [2(| and g aa = 2Huj±a aa as 
long as a aa < l±. 

Using Eq.0 for the coupling constants X q , we find 

a = *ma [— — V (10) 

where ^ a }, = mg a } 1 /h 2 p a . Similar to the discussion above, 
the interference between the direct interaction and the 
Rabi-transition terms in the Hamiltonian (0 leads to the 
disappearance of AQD-phonon coupling at = ir 2 /K 2 . 
The Luttinger liquid parameter K can thus be deter- 
mined by tuning the AQD to the decoupling point via 
a change in the known interaction constant j a f,. Thus, 
observation of the dynamics of the dot coupled to a Lut- 
tinger liquid provides a novel tool to measure K directly. 
For a weakly interacting ID liquid 7 aa <C 1 the Bogoli- 
ubov approximation applies, giving K = vr / ^/7 00 • In re- 
cent experiments 0] values 7 aa > 1 have been reached 
in in an array of independent one-dimensional tubes . 
In this case the critical value a c — 1 for the dissipative 
phase transition is reached at a a b/a a a ~ 1. 

In conclusion, we have shown that an AQD coupled 
to a superfluid reservoir leads to spin-Boson model with 
tunable parameters. The present model may be readily 
extended to arrays of AQDs, where the pseudo-spins can 
interact with "host" liquid in a collective manner. 
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